AD  7  3  4  7  9  8 


RADC-TR-71-249 
TECHNICAL  REPORT 
OCTOBER  1971 


EFFECTS  OF  TURBULENCE  INSTABILITIES 
ON  LASER  PROPAGATION 

RCA  LABORATORIES 


SPONSORED  BV 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
ARPA  ORDER  NO.  1279 


Approved  for  public  release; 
distribution  unlimited. 


The  views  end  conclusions  contained  in  this  document 
are  those  of  the  authors  and  should  not  be  interpreted 
as  necessarily  representing  the  official  policies,  either 
expressed  or  implied,  of  the  Advanced  Research  Projects 
Agency  or  the  U.S.  Government 


ROME  AIR  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
GRiFFISS  AIR  FORCE  BASE,  NEW  YORK 


When  U.S.  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than 
a  definitely  related  government  procurement  operation,  the  government  thereby  incurs  no  respon¬ 
sibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data  is  not  to  be  re¬ 
garded,  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  converying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  in  any  way  be  related  thereto. 


If  this  copy  is  not  needed,  return  to  RADC  (OCSE),  GAFB,  NY  13440. 


Security  Clanificitlon 


DOCUMENT  CONTROL  DATA  •  R  &  D 

(Sacurlty  elaaallteatlon  ol  tltla,  body  of  abaltael  and  Indaalng  annotation  muat  ba  antarad  whan  (ho  ovarall  rapott  la  elaaalllad) 

l  originating  activity  (Cotporata  author) 

RCA  Laboratories 

U.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

Princeton,  New  Jersey 

lb-°"OUP  N/A 

la.  REPORT  TITLE 

EFFECTS  OF  TURBULENCE  INSTABILITIES  ON  LASER  PROPAGATION 


4.  descriptive  NOTES  (Typ*  of  import  and  Inclualra  dalaa) 

_ Quarterly  Progress  Report  No.  1 _ 9  June  1971  -  8  September  1971 

a  authorisi  (FI  ft  namm,  mlddta  Initial,  laat  nama) 

David  A.  de  Wolf 


•  REPORT  date 

8  October  1971 


TO.  TOTAL  NO.  OP  PACES 

24 


7b.  NO.  OP  REPS 
11 


SO.  CONTRACT  OR  CRANT  NO. 

F30602-71-C-0356 

b.  PROJEC  T  NO. 

ARPA  1279 

e.  Program  Code  1E20 


SO.  ORIOINATOR'S  REPORT  NUMBER**) 


PRRL-71-CR-31 


9b.  OTHER  REPORT  NOlSl  (Any  other  numbers  dial  may  ho  aaatfnad 
I  th/e  roporfj 

j  RADC-TR-71-249 


10.  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

II.  SUPPLEMENTARY  NOTES 

12.  SPONSORING  MILITARY  ACTIVITY 

Monitored  by  Raymond  P.  Urtz,  Jr.; 

Advanced  Research  Projects  Agency 

(315)  330-3443  RADC  (OCSE)  GAFB, 

1400  Wilson  Blvd. 

NY  13440 

Arlington,  Va.  22209 

The  average  effective  area  of  the  focus  of  a  laser  beam  in  turbulent  air 
(diffraction-limited  in  free  space)  has  been  computed.  The  results  show  that 
the  focal  spot  can  be  decreased  only  to  a  certain  extent  by  increasing  the  trans¬ 
mitting  aperture  up  to  a  critical  size  determined  by  the  turbulence  structure 
constant  Cn2.  Critical-aperture  radius  roc  =  l/2kLl/2  (k  is  wavenumber, 

L  is  pathlength,  <m  is  inner-scale  wavenumber).  If  the  aperture  is  made  larger, 
the  focal-spot  area  will  fluctuate  around  a  constant  value  close  to  a  minimum 
value  L/kr0c  independent  of  the  aperture  size.  This  minimum  is  an  atmospheric 
limit  to  the  focussing  power  of  a  laser  for  applications  in  which  illumination 
is  not  shorter  than  a  typical  fluctuation  time. 


DD  .”".1473 


UNCLASSIFIED 

Security  Classification 


tfttCLftSsirifcp 

Security  Classification 


EFFECTS  OF  TURBULENCE  INSTABILITIES 
ON  LASER  PROPAGATION 

D.  A.  de  WOLF 


CONTRACTOR:  RCA  LABORATORIES 
CONTRACT  NUMBER:  F30602-71-C-0356 
EFFECTIVE  DATE  OF  CONTRACT:  9  JUNE  1971 
CONTRACT  EXPIRATION  DATE:  8  JUNE  1972 
AMOUNT  OF  CONTRACT:  $49,990.00 
PROGRAM  CODE  NUMBER:  1E20 

PRINCIPAL  INVESTIGATOR  AND 

PROJECT  ENGINEER:  DR.  DAVID  A.  de  WOLF 

PHONE:  609  452-2700  EXT.  3023 

CONTRACTING  ENGINEER: 

RAYMOND  P.  URTZ,  JR. 

PHONE:  315  330-3443 


Approved  for  public  release; 
distribution  unlimited. 


This  research  was  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department 
of  Defense  and  was  monitored  by  RADC 
Raymond  P.  Urtz,  Jr.  under  Contract  Number 


F30602-71-C-0356. 


! 

( 

PUBLICATION  REVIEW 

This  technical  report  has  been  reviewed  and  is  approved. 


FOREWORD 


This  Quarterly  report  was  prepared  by  RCA  Laboratories,  Princeton, 

New  Jersey  under  Contract  No.  F30602-71-C-0356.  It  describes  work  performed 
from  9  June  to  8  September  1971  in  the  Communications  Research  Laboratory, 

K.  H.  Powers,  Director.  The  principal  investigator  and  project  engineer  is 
D.  A.  de  Wolf. 

The  report  was  submitted  by  the  author  on  8  October  1971.  Submission  of 
this  report  does  not  constitute  Air  Force  approval  of  the  report's  findings 
or  conclusions.  It  is  submitted  only  for  the  exchange  and  stimulation  of 
ideas . 

The  Air  Force  Program  Monitor  is  Raymond  P.  Urtz,  Jr. 


iii 


ABSTRACT 


The  average  effective  area  of  the  focus  of  a  laser  beam  in  turbulent 
air  (diffraction-limited  in  free  space)  has  been  computed.  The  results  show 
that  the  focal  spot  cun  be  decreased  only  to  a  certain  extent  by  increasing 
the  transmitting  aperture  up  to  a  critical  size  determined  by  the  turbulence 
structure  constant  Cn  .  Critical-aperture  radius  roc  =  l/2kL^'2  Kffi1/6Cn 
(k  is  wavenumber,  L  is  pathlength,  xm  is  inner-scale  wavenumber).  If  the 
aperture  is  made  larger,  the  focal-spot  area  will  fluctuate  around  a  constant 
value  close  to  a  minimum  value  L/kroc  independent  of  the  aperture  size.  This 
minimum  is  an  atmospheric  limit  to  the  focussing  power  of  a  laser  for  applica¬ 
tions  in  which  illumination  is  not  shorter  than  a  typical  fluctuation  time. 
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SECTION  I 


SUMMARY 


When  images  are  formed  from  laser  beams  propagating  through  turbulent 
air,  a  variety  of  scintillation  phenomena  or  curs:  beam  wander,  intensity 
fluctuations,  hot-  and  cold-spot  formation,  image  blurring,  spot  broadening, 
etc.  The  purpose  of  this  project  is  to  study  these  effects  analytically, 
and  thus  to  interpret  measurements  and  predict  performance  in  future  laser 
sys  terns .  (_ 

The  difficulties  of  an  analytic  approach  -  aside  from  obvious  mathemati¬ 
cal  ones  -  lie  in  the  present  status  of  theoretical  knowledge  of  the  basic 
optics  in  turbulent  air.  The  statistics  of  the  irradiance  (the  prime  quantity 
of  interest)  are  uncertain,  and  in  particular  its  variance  and  higher-order 
measures  of  fluctuations  are  not  clearly  established. 

Rather  than  starting  a  basic  investigation  in  this  first  quarter  into 
the  full  statistics  of  the  irradiance  of  an  Idealized  wave,  we  preferred  to 
define  and  solve  a  more  limited  problem:  the  limitations  imposed  by  atmo¬ 
spheric  turbulence  upon  the  focussing  of  a  laser  beam.  It  is  quite  clear 
that  a  primary  factor  in  reconnaissance  and  weapons  performance  of  laser 
beams  is  the  ability  to  focus  the  beam.  In  contrast  to  the  non-linear  and 
self-induced  effects  such  as  thermal  blooming  and  distortion,  atmospheric 
turbulence  creates  beam  distortionc  that  are  not  a  function  of  beam  strength. 
Furthermore  the  effects  are  random  and  therefore  unpredictable. 

Previous  work  on  this  problem  appeared  to  us  to  be  incomplete,  or 
-  perhaps  -  not  correct.  We  have  defined  an  effective  focal-spot  area  of 
a  laser  beam,  and  have  computed  its  average  area  in  the  presence  of  atmospheric 
turbulence.  The  analysis  utilizes  only  those  irradiance  statistics  that  are 
well-established  and  agreed  upon  by  most  workers  in  this  field.  The  results 
indicate  that  for  sufficiently  small  aperture  radii  rc  one  can  get  large 
focal  spots  with  radius  r^  ^  L/krQ  (L  is  pathlength,  k  is  wavenumber)  which 
can  be  made  smaller  by  increasing  rQ  until  the  aperture  radius  approaches  a 
critical  size  determined  by  the  atmospheric-turbulence  strength  parameter 
Cn^,  namely  roc  =  l/2kLl(/2|<ml/6cn  (k  is  wavenumber,  L  is  pathlength,  is 
inner-scale  wavenumber) .  There  may  be  no  advantage  in  increasing  r0  beyond 
this  critical  size  because  the  focal  spot  then  fluctuates  in  size  and  on  the 
average  does  not  decrease  in  radius  below  L/kroc.  Details  of  the  analysis 
are  appended  in  technical-report  form  in  subsequent  sections  of  this  Quarterly 
Report. 

The  results  -  given  for  so-called  Kolmogorov  turbulence  -  can  be  tested. 
Comparisons  with  data  -  to  be  made  available  yet  -  are  contemplated  in  the 
future.  If  necessary,  the  results  can  be  adapted  to  other  forms  of  turbulence. 
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SECTION  II 


FOCUSSED  LASER  BEAM  IN  TURBULENT  AIR:  INTRODUCTION 


The  topic  of  this  report  is  the  effect  of  turbulent  air  upon  the  focussing 
of  a  laser  beam  aimed  horizontally  within  a  few  meters  of  earth's  or  sea's 
surface  of  a  receiver  located  a  distance  L  from  the  effective  optical  aperture 
(with  radius  r^)  of  the  beam-producing  laser  system. 

In  free  space,  the  focal  spot  is  diffraction  limited.  Its  shape  is 
dependent  upon  the  way  the  illuminating  flux  varies  over  the  aperture  located 
in  the  plane  z  =  0.  Common  to  all  practical  cases  is  the  property  that  the 
focal  spot  has  an  effective  radius  determined  by  the  distance  r^  ■  L/kr0 
(k  is  the  wavenumber  of  the  laser-beam  radiation,  which  may  be  considered 
monochromatic  without  undue  restriction).  It  cannot  be  made  smaller  without 
modifying  L,  k,  or  rfl.  In  turbulent  air,  however,  the  phase  relationships 
between  rays  emanating  from  different  parts  of  the  transmitting  aperture  to 
the  focus  are  modified  and  the  result  is  spot  broadening  beyond  the  diffrac¬ 
tion  limit. 

A  naive  calculation  of  the  effect  can  be  given  first.  Consider  the 
aperture  of  the  laser-optics  system  producing  a  focussed  beam  as  a  lens. 

By  exploiting  the  refractive  properties  of  the  lens,  i.e.,  by  choosing  its 
material  and  shape,  one  attempts  to  direct  all  the  rays  penetrating  it  to  a 
focal  point.  However,  ray  bending,  or  diffraction  of  rays,  occurs  around  the 
edges  of  the  lens  and  deflection  of  rays  by  an  angle  8  ^  (kr0)-1  occurs.  As 
a  result,  light  from  the  edges  of  the  lens  which  is  intended  for  the  focal 
point  on  the  central  axis  at  distance  L  is  misaimed  by  an  angle  of  the  order 
of  0.  A  blurring  of  the  focal  point  into  a  spot  of  thickness  rj,  ^  L0  occurs. 

Of  course  this  picture  is  only  one  of  a  number  of  other,  equivalent,  ways  of 
describing  diffraction.  Now,  when  the  rays  must  propagate  through  a  random 
geometrical-optical  atmosphere  (i.e.,  a  medium  which  is  quite  transparent  but 
which  deflects  rays  continuously  with  local  radius  of  curvature  very  large 
compared  to  wavelength  in  randomly  varying  directions)  the  misaiming  is  by  an 
angle  0  +  60,  and  the  blurring  occurs  over  a  thickness  rj,  L(9  +  60).  Be¬ 
cause  the  average <L  60>  =  0,  the  average  blurring  radius  in  one  direction  is 
as  before:  diffraction  limited.  The  average  area  ir  <rj^>  is  a  measure  of 
the  net  effect  of  the  atmosphere,  and  it  is  easily  observed  that 

<rL2>  =  L2  [o2  +  <(60) 2>]  (1) 

2 

The  angular  variance  <(60)  >  in  homogeneously  turbulent  air  can  be  estimated 
by  geometrical-optics  formulas  for  ray  curvature.  Although  the  calculation 
is  not  new,  it  is  not  easily  accessible  and  we  provide  a  short  derivation  in 
Appendix  II.  We  also  set  0  =  (krQ)“l  and  find 

< r  2>  *  (L/kr  )2  [l  +  3.9  C  2  k2  L  k  1/3  r  *1  (2) 

L  o  L  n  m  o  j 


2 


This  formula  predicts  that  the  effects  of  turbulence  are  negligible  unless 
Cn^  approaches  a  critical  value  determined  by  the  parameters 


K  -1'3 


(3) 


In  words,  the  broadening  effect  of  turbulence  is  increased  by  utilizing 
higher  frequencies  and  larger  apertures.  Lawrence  and  Strohbehn[l] ,  however, 
have  utilized  another  type  of  naive  argument  to  predict  a  critical  Cn  ,  and 
their  result  amounts  to 


crit. 


(4) 


In  contrast  to  the  first  line  of  reasoning,  this  result  predicts  a  frequency- 
independent  broadening  that  decreases  with  increasing  aperture.  There  are 
perhaps  some  reasons  to  be  apprehensive,  especially  of  this  latter  line  of 
reasoning,  but  a  simplified  argument  such  as  the  above  has  little  value  other 
than  to  provide  an  intuitive  interpretation  for  a  rigorous  derivation.  Tt 
should  be  accompanied  by  a  rigorous  derivation.  In  the  following  sections, 
Equation  (2)  will  be  derived  rigorously  for  a  focussed  laser  beam  in  turbulent 
air  governed  by  a  Kolmogorov  spectrum  of  refractive-index  fluctuations  for 
propagation  paths  that  exceed  those  for  which  the  irradiance  saturates [2 ] . 


3 


SECTION  III 

FOCUSSED  LASER  BEAM  IN  TURBULENT  AIR:  BASIC  FORMULATION 


Consider  the  aperture  of  a  laser  system  in  the  z  =  0  plane.  A  point  on 
the  aperture  is  denoted  by  r^  =  (p"i,0).  Schmeltzer [3]  has  shown  that  the  elec¬ 
tric  field  EQ(r)  at  location  r  =  (p",L)  in  free  space  is  given  by  the  aperture 
integral 

Vj)  ’  -  fc  if  Jd2pl  Go  (™l)  Uo  (*l)  •  (5) 

Here,  GQ(Ar)  =  [exp(ikAr)]  /Ar,  is  the  field  of  a  spherical  wave  emanated  at 
an  origin  and  received  at  a  vectorial  location  tr  =  r-r^  from  that  origin;  R 
is  the  radius  of  curvature  of  the  transmitted  wavefront,  and  UQ  is  the  aperture 
function  weighting  the  amplitude  of  the  transmitted  field  [it  is  chosen  as  a 
Gaussian  in  Equation  (5)  to  conform  with  laser  beams,  but  it  can  also  be 
chosen  differently  to  correspond  to  other  optical  systems  that  produce  a 
positively  or  negatively  (divergent)  focussed  beam] .  It  can  be  seen  that  the 
approximation  9GQ/3z  %  ikGQ  is  excellent  in  Equation  (5)  .  This  approximation 
ignores  terms  of  order  (kL)“l  and  of  order  L/krQ2  in  the  electric  field;  it 
does  not  ignore  the  phase  effects  of  the  exponent  of  Gq. 

The  formulation  [Equation  (5)]  with  9G  /9z  replaced  by  ikGQ  adapts  it¬ 
self  excellently  to  the  problem  of  a  laser  beam  in  turbulent  air.  The  fact 
that  Equation  (5)  is  a  linear  integral  equation  with  kernel  GQ  over  an  effec¬ 
tive  source  function  immediately  establishes  that  in  turbulent  air 

E(r)  -  -  §■  J d2px  G  (r.tjJ  Uq  (rj)  exp^-ikp^m)  ,  (6) 

where  G(r,r^)  is  the  spherical-wave  field  at  r  when  a  spherical  wave  is 
emitted  from  rQ  in  turbulent  air.  The  field  G(r,r^)  reduces  to  GQ(r-ri)  as 
the  turbulence  vanishes  and  the  medium  is  thus  reduced  to  free  space,  conse¬ 
quently  we  find  it  useful  to  set  G(r,r^)  h  B(r ,r^)G0(r-r^)  so  that  B  -►  1 
as  6e  **■  0. 

If  the  beam  is  focussed  at  L,  i.e.,  R  =  L  one  can  approximate  G^r-r^) 
to  first  order  in  the  exponent  and  to  zero-th  order  in  the  denominator  to 
obtain, 

E(r)  =  L-1  exp[ik(L  +  p 2/2l)]  •  ^  ,rl)Uo (rl)  exp(ikp -P]/1)  •  (7) 
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This  formulation,  which  also  ignores  terms  of  order  kr0^/L^ 
for  our  purnose.  The  received  power  at  r  is  proportional 
form  the  ir radiance  I  =  | E | ^  to  find, 


is  very  convenient 
| E | 2 .  We  therefore 


I(r)  =  |  (k/2irL)  J  d^  B^r.r^  U^r^  exp  (ikp  .p^L^  2  (8) 


The  irradiance  distribution  over  the  focus,  i.e.,  the  dependence  of  I(r)  upon 
p  for  fixed  L,  determines  the  properties  of  the  focal  spot.  The  free-space 
properties  are  immediately  apparent  from  Equation  (8)  because  B  =  1.  In  that 
case,  the  insertion  of  Equation  (5)  for  U0(ri)  yields  an  easily  performed 
Fourier  integral  inside  the  modulus  signs  of  Equation  (8) : 


10M  - 


(9) 


and  it  can  be  seen  by  making  the  integration-variable  p^,  non-dimensional  for 
more  general  UD(r^)  that  IQ(r)  depends  only  upon  krQp/L,  and  so  that  IQ(r) 
decreases  (albeit  not  monotonically)  with  increasing  krQp/L.  It  is  obvious 
from  Equation  (9)  and  its  generalization  that  the  focal  spot  has  a  halfwidth 
rL  that  is  easily  related  to  L/krQ.  The  difficulty  is  in  extending  the  con¬ 
cept  to  the  case  Bj*l,  i.e.,  to  turbulent  air,  because  it  is  not  possible  in 
general  to  perform  the  random  Fourier  integral. 
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SECTION  IV 


DEFINITION  OF  FOCAL-SPOT  RADIUS  IN  TURBULENT  AIR 


As  with  so  many  other  randomly  varying  quantities,  one  can  only  hope  to 
extract  more  information  about  the  behavior  of  I(r)  with  the  atmospheric  and 
system  parameters  by  calculating  statistics.  The  random  quantity  I(r)  In 
Equation  (8)  is  determined  by  the  statistics  of  B(r ,?i)B*(r,?2)  at  two  aper¬ 
ture  locations  r^  and  r2 •  The  statistics  of  the  normalized  spherical-wave 
field^B(r ,ri)  are  uncertain  to  date,  but  the  average  of  the  above  product, 
<B(f,ri)B*(r,r2)>  the  mutual  coherence  factor  (mcf),  has  been  determined  by 
a  variety  of  methods [4-6]  for  long  propagation  distances. 

°° 

<B(r,r1)B*^r,r2^>  -  exp 


Ap  =>  Ipj-f^j 

Here,  JQ(KAp)  is  a  Bessel  function  (zero-th  order),  and  the  reciprocity 
theorem  has  been  invoked  to  interpret  B(r,r^)  also  as  the  normalized  field 
for  propagation  from  ?  to  r^.  The  validity  of  Equation  (10)  is  limited  only 
by  the  conditions  e  2  <<  fc-2  j^-2  antj  ^  <<  j^e”  aside  from  the  usual  optical 
(far-field  and  sagittal)  assumptions.  The  region  of  validity  includes  that 
in  which<  [I(?)  ]2>  is  observed  to  saturate. 

Thus,  <I(r)>  appears  promising  as  a  statistic  which  yields  information 
about  rL.  Unfortunately  the  computation  of  <I(r)>  after  inserting  Equation 
(10)  into  the  average  of  Equation  (8)  is  still  very  difficult.  Consider 
Equation  (10) .  For  propagation  over  more  than  several  meters  under  at  least 
moderately  turbulent  conditions  (Cn2  >  10“^^m”^'^)  the  quantity  k2  e2L  1. 
It  therefore  follows  that  the  coordinate  Ip-^— p£ | >  which  we  abbreviate  by  Ap, 
need  be  specified  in  the  integrand  only  for  Ap  <  L0  because  <B(1)B*(2)> 
(another  obvious  abbreviation)  becomes  negligibly  small  otherwise.  For 

<Ap<  Lq  it  then  follows  that  the  exponent  in  Equation  (10)  is  proportional 
to  k2c2L(Ap)-*/3 ,  and  to  k2Cn  Licml/3(Ap)2  for  ap  <  £q.  Even  so,  it  appears 
that  <I(r)>  can  be  evaluated  further  only  numerically. 

However,  the  form  of  Equation  (8)  suggests  strongly  that  an  integration 
over  d2p,  the  plane  of  the  focal  spot,  will  simplify  matters  because  a  delta 
function  in  Ap  (or  something  like  it)  follows  when  B(1)B*(2)  depends  weakly 
upon  r.  The  first  radius-like  quantity  we  can  form  is 

<rL2>  i  f  d2p  <1  (r)>  /  < I  (0,L)>  (11) 


1  -#1/  dKK*<K>  C1  -  L_1f 


dz  JQ(2KzAp/L) 


II 


(10) 
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The  numerator  is  easily  evaluated  -  we  shall  do  so  presently  -  but  the  de¬ 
nominator  is  subject  to  the  problems  discussed  in  the  preceding  paragraph. 
Hence  we  define  an  effective  radius  with  the  next  non-zero  higher  order 
weighted  mean: 


> 


f 


d2p  p2  <I(p,L)> /f i i 


d2p<I(p,L)> 


(12) 


The  denominator  of  Equation  (12)  is  easily  computed  because  the  d  p  integra¬ 
tion  introduces  a  two-dimensional  delta  function  in  ^“^2*  Consequently 
<B(1)B*(2)>  reduces  to  <B(1)B*(1»  =  1.  What  is  left  yields 


2 

just  as  we  expect  from  energy  conservation  (irr0 
the  focus  by  the  aperture) . 


(13) 


is  the  power  radiated  toward 


Tne  numerator  of  Equation  (12)  requires  somewhat  more  work.  After  per¬ 
formance  of  the  d^p  integral,  we  obtain 


<B(1)B*(2)> 


exp 


[-  (°12  +  c22V2ro2] 


x 


+ 


(14) 


We  perform  partial  integration  with  respect  to  x^  and  y^  twice  in  order  to 
remove  derivatives  from  the  delta  function.  Because  the  boundary  terms  at 
xi  =  +  °°,  yi  =  +  00  ate  zero,  the  net  effect  is  to  transfer  the  differential 
operator  in  Equation  (14)  to  the  two  other  factors  of  the  integrand.  The 
presence  of  a  delta  function  then  requires  that  we  know<B(l)B*(2)>  only  for 
infinitesimal  (but  non-zero)  Ap  in  order  to  obtain  the  derivatives  and  then 
set  Ap  *  0.  This  is  a  considerable  simplification.  Referring  to  Appendix  I, 
we  observe  that 


lim  <  B(1)B*  (2)> 

-  lim 

exp 

Ap+0 

Ap+0 

*  lim 

exp 

Ap-*0 

^(k2e2L/24ir)  M3.  (aP)2  ] 

[-  A  ‘a1'3  C„2  W2] 


(15) 
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Here,  we  have  utilized  Equations  (1-6)  and  1-7)  in  the  second  step.  Let  us 
define  B  =  1.3  k^LKm^' ^Cn^,  so  that  Equation  (14)  reduces  to 


(16) 


The  differentiations  are  now  easily  performed,  and  after  then  setting  ip  *  0 
we  find  j 


2 

P  P 


2 


-  ir(L/k)2  (l  +  4BrQ2  ) 


(17) 


The  averaged  effective  rrlius  defined  in  Equation  (12)  is  thus  given  by 
Equations  (13)  and  (17),  and  the  result  is 


<rL2>  =(L/kro)2  (X  +  48ro2) 


6  5  1.3 


k2  L  k  1/3 
m 


(18) 


It  is  interesting  to  note  that  this  is  roughly  equal  to  Equation  (2)  [the 
numerical  coefficient  of  k^Lic^^C  * (Ap) *  is  3  y  (*)  y  (e2)/3i  3.9  in  the 

former  case] .  The  result  shows  the  following:  Given  a  laser  system  (L  and  k 
specified)  and  atmospheric  conditions  (Cn^  specified),  it  follows  that  the 
focal  spot  can  be  made  smaller  by  increasing  rQ  until  it  approaches  a  critical 
value  roc  =  1/2B  There  is,  however,  no  further  decrease  of  average 

focal-spot  area  when  increasing  rQ  above  roc,  because  <r^>  saturates  rapidly 
to  the  value  (L/kroc)^  which  is  independent  of  ro« 

The  derivation  can  be  extended  to  non-Gaussian  U(?^)  by  defining  an 
aperature  radius 


}o/S*2>\ 


U 

o 


(19) 


The  same  result  is  obtained  as  before,  namely  Equation  (18)  with  the  above 
new  definition  of  rQ,  and  the  derivation  does  not  differ  essentially  from 
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that  for  Gaussian  UD(r^) .  The  only  requirement  is  that  r  as  defined  in  Equation 
(19)  is  finite.* 
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SECTION  V 


CONCLUSIONS 


In  the  previous  section  we  derived  rigorously  an  expression  for  an  irra- 
diance-weighted  focal-spot  area  on  the  average.  This  area  is  broadened  by  at¬ 
mospheric  turbulence  beyond  the  diffraction-limited  area  ir  (L/krQ)2.  Surpris¬ 
ingly  enough,  the  simple  argument  of  Section  I  yielded  exactly  the  same  answer 
although  it  ignores  correlation  of  rays  from  different  parts  of  the  transmit¬ 
ting  aperture  and  utilizes  the  ray  concept  as  if  there  were  no  ray  crossings 
and/or  caustic  effects  before  the  focal  point  is  reached.  We  would  like  to 
understand  this  somewhat  more.  Consider  the  electric  field  E(p,L)  received  in 
the^focal  plane,  and  Fourier  analyze  it  in  terms  of  transverse  random  components 
dE(K,L) 

E(p,L)  *=  / dE(K,L)  exp(-iK.p) 

with  <dE(K1,L)dE*(K2,L)>  *  (4ir)~1d2K1d2K2  (20) 


This  is  tjje  customary  Fourier  analysis  for  a  homogeneous  isotropic  random  pro¬ 
cess  and  I(K,L)  is  the  two-dimensonal  spectral  density  of  E(p,L).  From  Equation 
(20)  and  the  definition  (12)  it  is  then  not  difficult  to  find 


<rL2> 


-  [a^(«.L)/«2]K.0/ 1  (o.L) 


(21) 


i.e.,<rLz>  is  determined  by  smallest-wavenumber  components  of  the  spectral 
density  of  E(p,L)  in  the  focal  plane.  The  main  effect  determined  by  small- 
wavenumber  m£,L)  can  only  be  due  to  a  gross  deflection  of  the  entire  beam. 
Apparently  this  is  why  the  effect  of  atmospheric  turbulence  on  <xi^>  is  so 
aptly  described  by  a  ray  angle  60  as  in  Section  I. 

If,  then,  the  quantity  <r^^>  is  any  measure  of  the  average  effective  area 
of  the  focal  spot  in  turbulent  air,  one  finds  that  the  turbulence  strength  Cn2 
must  exceed  a  critical  value  given  parametrically  in  Equation  (3)  in  order  for 
appreciable  broadening  beyond  the  diffraction  limit.  This  differs  from  the 
conclusion  of  Lawrence  and  Strohbehn[l]  in  Equation  (4).  A  corollary  of  our 
result  is  that  the  only  restriction  on  aperture  size  for  validity  of  the  approx¬ 
imations  is  L  <<  krQ2.  We  have  used  this  restriction  in  developing  Equation 
(6) .  It  has  also  been  used  in  ignoring  the  fact  in  Equation  (10)  that  r  makes 
an  angle  0r  £  r/L  ^  (kr0)-^  <<  (kL)“^/‘  with  the  central  axis  (terms  of  order 
k_^L~l  are  discarded).  The  results  [Equation  (18)]  are  displayed  in  Figures  1 
and  2  for  wavelengths  X  =  0.6  pm  and  X  =  10.6  um  in  the  following  fashion:  We 
defined  a  critical-aperture  radius  roc  =  (40) 2  so  that  for  given  Cn2  and 
L  the  average  focal-spot  area  is  just  twice  the  diffraction  area.  Let 
r^  'v  L/kroc.  If  rQ  <  roc  then  the  diffraction- limited  radius  ^  L/krQ  is 
a  good  approximation,  and  rL  >  rLc.  When  rQ  >  roc,  then  Equation  (18)  shows 
that  rL  is  slightly  larger  than  rLc.  Consequently  <r^  :->is  the  smallest 
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L  in  km  — ► 


Figure  1.  Maximal  aperture  radius  roc  required  to  obtain 
a  minimal  focal-spot  area  with  effective  radius 
rLC.  Radii  roc  and  rLC  as  functions  of  path- 
length  L  for  X  =  0.6  ym,  <m-^  =  7  mm.  and 
Cn2  =  10_14(MT),  10_15(ST),  and  10_l6(VST)  in 
m-  2/3  units. 


average  area  of  the  focal  spot  to  be  obtained  under  atmospheric  conditions 
prescribed  by  given  values  of  C  .  The  definitions  yield 


r  -  l/2kL1/2<  1/6  C 
oc  m  n 


Lc 


2L3/2.  1/6  C 
m  n 


(22) 


and  these  quantities  are  plotted  for  <m  %  7  mm  in  Figures  1  and  2  for  three 
values  of  C  2  corresponding  to  moderate  (Ml),  strong  (ST),  and  very  strong 
(VST)  turbulence. 


Finally,  we  compare  our  results  to  previous  work.  Gebhardt  and 
have  computed  rT /r  vs  kr  ^/L.  Their  results  show  the  same  trend  as 

L  O  O 


Collins[7] 
ours  for 
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L  in  km - ► 

Figure  2.  Maximal  aperture  radius  roc  required  to  obtain 
a  minimal  focal-spot  area  with  effective  radius 
rLc.  Radii  roc  and  rLc  as  functions  of  path- 
length  L  for  A  =  10.6  ym,  Km“l  *  7  mm,  and 
Cn2  =  10-J4(MT),  10_15(ST),  and  l<r16(VST)  in 
m“2/3  units. 


rQ  «  rQC.  but  there  is  no  saturation  for  larger  rQ.  As  a  consequence,  they 
conclude  -  in  sharp  contrast  to  us  -  that  better  focussing  can  be  achieved  by 
increasing  rQ  no  matter  how  large  Cn^  is.  The  difference  in  their  work  appears 
to  arise  from  the  use  of  weak-amplitude  theory  and  a  number  of  assumptions  re¬ 
gard  *.g  the  average  log-amplitude.  Lutomirski  and  Yura [8]  compute  a  half-power 
anglv-  ^  l.l/2Kml/6cn  *s  similar  to  our  result  for  ri/L  (also  a  half-power 

angle)  when  r„  >  roc.  The  result  appears  to  agree  with  ours  only  for  large 
paths  L  »  k”2cn_2Km_^' 3 f  however.  Their  choice  of  aperture  (equivalent  to 
that  discussed  in  the  footnote  on  page  9)  may  make  further  comparison 
difficult. 

Note  added  in  proof:  after  the  completion  of  this  work,  it  came  to  our 
attention  that  Varvatsis  and  Sancer[ll]  have  found  the  same  result. 
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APPENDIX  I 


TURBULENCE  QUANTITIES  REQUIRED  IN  WAVE-PROPAGATION  CALCULATIONS 


Different  authors  use  differing  nomenclature  and  normalizations  of  turbu¬ 
lence  scales  and  spectra.  To  avoid  confusion,  we  tabulate  our  own  definitions 
for  future  reference. 

The  basic  quantity  is  the  deviation  from  the  mean  dielectric  permittivity 
6e(r).  Its  lowest  non-zero  statistic  is  its  variance  [6e(r)]^  which  may  or 
may  not  be  dependent  upon  location  r.  In  the  case  that  it  is  not,  we  abbreviate 
the  variance  by  the  symbol  ^ . 

In  the  case  of  homogeneous  turbulence,  the  twc-point  correlation 
<6c(rj)  6c  (r2)>  depends  only  upon  the  difference  vector  Ar  =  r^  ~*2‘ 

We  write  this  two-point  correlation  as  ^C(Ar),^thus  insuring  the  normalization 
C(0)  =  1,  and  define  the  turbulence  spectrum  $(K)  as 


<J(K)  =  f  d3Ar  C  (Ar)  exp(iK.Ar)  (1-1) 


When  there  is  no  preferred  direction^(isotropic  turbulence) ,  the  spectrum 
depends  only  on  the  length  of  the  vector  K,  and  we  will  denote  the  argument 
as  a  scalar  quantity.  For  some  purposes,  we  prefer  to  give  special  attention 
to  the  direction  of  propagation:  the  z  axis  in  our  work.  Vector  K,  for  example, 
will  then  be  denoted  as  K^,  Kz  in  functional  arguments  (where  KT  is  the  projec¬ 
tion  of  K  upon  the  z  =  0  plane) . 

Definition  (1-1)  implies  for  the  case  of  isotropy  that 
00 

M2  E  j  dKK2$(K)  -  2tt2  (1-2) 

0 

The  Kolmogorov  spectrum  $(K)  corresponding  to  the  "two-thirds"  law  of  the 
temperature  structure  function  in  the  inertial  subrange  of  turbulence,  and 
modified  to  incorporate  small  and  large  wavenumber  properties  is 

$(K)  =  y($)  Lq3  (l  +  K2Lo2)~11/6  exp(-K2/Km2)  (1-3) 

L  :  macroscale  of  turbulence 

£°:  microscale  of  turbulence  (£Q  =  5.92  icm  ) 

Y(<J?:  a  normalizing  constant  obtained  from  Equation  (A-2)  and  given  to  good 
accuracy  by  47r2r(ll/6) /T(3/2)  1(1/3)  *  15.7. 
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An  important  parameter  is  the  integral  scale  defined  as  the  first  moment 
M-^  divided  by  16  n 

00 

Vlfe  /  dKK*(K>  -*(li)Lo  (I‘4) 

0 

a  constant  obtained  by  inserting  Equations  (1-3)  into  (1-4)  and  given  to 
good  accuracy  by  3  y($)/80ir  «  1.88 

For  some  applications,  it  will  be  useful  to  have  an  expression  for  the  first 
moment  itself : 


00 


0 

2  2 

There  is  a  connection  between  e  and  the  structure  constant  Cn  .  By  forming 
the  refractive-index  structure  function  D(r)  by  Tatarski's[9]  Equations  (1.38) 
we  obtain  in  terms  of  our  normalizations 


D(r) 


16n 


J d3K  $(K)  [l-cos(ic.r)  J 


with  <KK)  given  by  Equation  (1-3).  We  compare  this  to  Tatarski 's [9]  Equations 
(1.41)  and  (3.24)  and  obtain 


*2  5  y  (,2)  C  2L  2/3 
'  '  n  o 

Y  (e2)  =  32tt3  x  0.033/y($)  *  2.08 


(1-7) 
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APPENDIX  II 


RAY  BENDING  IN  TURBULENT  AIR:  ANGULAR  VARIANCE 

2 

A  brief  derivation  of  the  angular  variance  < (66)  >  for  ray  bending  over 
distance  L  in  turbulent  air  is  given  here.  In  geometrical  optics [11]  the 
local  radius  of  curvature  is  tlje  inverse  length  of  the  vector  n“l  ^n  where 
n  is  the  refractive  index  and  Vj  is  the  gradient  in  the  plane  through  the 
local  portion  of  the  ray,  normal  to  the  local  propagation  direction.  In 
turbulent  air,  this  vector  is  approximated  well  by  V^.6e/2,  because 
n  =  (e)l'2  «  i  +  <5e/2  with  6e  «  1.  The  geometry  of  Figure  II-l  shows  that 
ds  =  Rcd60  and  consequently 

L  L 

«0  -  f  ds  rJrc(s)  =  |  I  dsVT6e(s)  (II-l 

0 


>s 

The  vector  notation  (Rc  is  a  unit  vector)  indicates  an  integration  over  direc¬ 
tion  of  60  (i.e.,  over  change  of  plane  in  which  66  lies)  as  well  as  over  the 
change  in  curvature.  We  form  <(60)  from  (II-l) ,  replace  the  gradient  opera¬ 
tors  at  points  s^  and  S£  of  the  double  integral  by  the  gradient  operator  with 
respect  to  s^-S2«  then  apply  the  formulas  of  Appendix  I  to  find 

QO 

<  ( 6  6 ) 2  >  =  (e2L  /  8tt)  j  dKK34>(K) 

0  (II-2) 

«  f*3y  (5>) y  (e2)/ 8tt~]  C  \  k  1/3  «  3.9  C  k  1/3 

L  'Jn  m  n  m 

The  above  equation  is  a  geometrical-optics  estimate  of  the  mean  square  angular 
deviation  of  a  ray  in  turbulent  air  with  a  modified  Kolmogorov  spectrum  of 
refractive- index  fluctuations. 
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